Scale-free Networks on Lattices 



O 

o 

(N 



I 

C/3 



C3 



I 

C 

O 

o 



> 

m 

in 
o 

<N 
O 

c3 



c 

o 
o 



X 
S3 



Alejandro F. Rozenfeld 1 , Reuven Cohen 1 , Daniel ben-Avraham 2 , and Shlomo Havlin 1 

1 Minerva Center and Department of Physics, Bar-Ran University, Ramat-Gan 52900, Israel 
2 Department of Physics, Clarkson University, Potsdam, New York 13699-5820 

(February 1, 2008) 

We suggest a method for embedding scale-free networks, with degree distribution P(k) ~ k~ x , 
in regular Euclidean lattices. The embedding is driven by a natural constraint of minimization of 
the total length of the links in the system. We find that all networks with A > 2 can be successfully 
embedded up to an (Euclidean) distance £ which can be made as large as desired upon the changing 
of an external parameter. Clusters of successive chemical shells are found to be compact (the fractal 
dimension is df = d), while the dimension of the shortest path between any two sites is smaller than 
one: d m 



X-l-l/d 



contrary to all other known examples of fractals and disordered lattices. 



Many social, biological, and communication systems 
can be properly described by complex networks whose 
nodes represent individuals or organizations and links 
mimic the interactions among them Q. An important 
class of complex networks are the scale-free networks, 
which exhibit a power-law connectivity distribution. Ex- 
amples of scale-free networks include the Internet [ 
WWW @||, metabolic |j| and cellular networks 
Most of the work done on scale free networks concerns 
off-lattice systems (graphs) where the Euclidean distance 
between nodes is irrelevant. However, real-life networks 
are often embedded in Euclidean space (e.g., the Internet 
is embedded in the two-dimensional network of routers, 
neuronal networks are embedded in a three-dimensional 
brain, etc.). Indeed, in the case of the Internet, indi- 
cations for the relevance of embedding space is given in 



In this Letter we develop a method for generating scale- 
free networks on Euclidean lattices and study some of its 
properties. As a guiding principle we impose the natural 
restriction that the total length of links in the system be 
minimal. 

Our model is defined as follows. To each site of a d- 
dimensional lattice, of size R, and with periodic bound- 
ary conditions, we assign a random connectivity k taken 
from the scale-free distribution 



P(k) = Ck 



-x 



<k<K, 



(1) 



where the normalization constant C « (A — l)m A_1 (for 
K large) || . We then select a site at random and connect 
it to its closest neighbors until its (previously assigned) 
connectivity k is realized, or until all sites up to a dis- 
tance 



r(Jfc) = Ak 1/d 



(2) 



have been explored. (Links to some of the neighboring 
sites might prove impossible, in case that the connectivity 
quota of the target site is already filled.) This process is 
repeated for all sites of the lattice. We show that follow- 
ing this method networks with A > 2 can be successfully 
embedded up to an (Euclidean) distance £ which can be 



made as large as desired upon the changing of the exter- 
nal parameter A. 

Suppose that one attempts to embed a scale-free net- 
work, by the above recipe, in an infinite lattice, R — > oo. 
Sites with a connectivity larger than a certain cutoff 
k c (A) cannot be realized, because of saturation of the 
surrounding sites. Consider the number of links n(r) en- 
tering a generic site from a surrounding neighborhood of 
radius r. Sites at distance r' are linked to the origin with 
probability P(k' > (r'/A) d ): 



P\k' > 




= C 



k~ A dk 



(3) 



Hence 



n(r) 



(5r) d 

1 r' < A. 

(^)d(l-A) r i > A 



dr'r' d - x P k' > 




(4) 
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d{\ - 2) d(X - 2] 
The cutoff connectivity is then 

1 



k c = lim n(r) 



A 



.4" 



The cutoff connectivity implies a cutoff length 
£ - r(k c ) ~ (A - 2)- 1 ' d A 2 . 



(5a) 



(5b) 



The embedded network is scale-free up to distances r < £, 
and repeats itself (statistically) for r > £, similar to the 
infinite percolation cluster above criticality: The infinit 
cluster in percolation is fractal up to the coherence length 
£ and repeats thereafter Jl0|~p^[ . 

When the lattice is finite, R < oo, the number of sites 
is finite, N ~ R d , which imposes a maximum connectiv- 

ity @01 



l 



This implies a finite-size cutoff length 
r moa = r(J0~AiZ 1 /< A - 1 >. 



(6a) the distribution terminates at the cutoff k c . The scale- 
free distribution is altered slightly, for k < k c , due to sat- 
uration effects, but the overall trend is highly consistent 
with the original power-law. The scaling in the inset con- 

(6b) firms that k c ~ A d . In Fig. 0(b), £ > r max and the cutoff 



The interplay between the three length scales, R, £, r max , 
determines the nature of the network. If the lattice 
is finite, then the maximal connectivity is k max = K 
only if r max < £. Otherwise (r max > £) the lattice re- 
peats itself at length scales larger than £. As long as 
min(r max ,£) <C R, the finite size of the lattice imposes 
no serious restrictions. Otherwise ( m.m(r max , £) > R) 
finite-size effects become important. We emphasize that 
in all cases the degree distribution (up to the cutoff) is 
scale- free. 

In Fig. |](a) we show typical networks that result 
from our embedding method, for A = 2.5 and 5 in two- 
dimensional lattices (in this Letter, we limit our numeri- 
cal results to d — 2). The larger A is the more closely the 
network resembles the embedding lattice, because longer 
links are rare In Fig. [l](b) we show the same net- 
works as in part (a) where successive chemical shells are 
depicted in different shades. Chemical shell I consists of 
all sites at minimal distance (minimal number of con- 
necting links) / from a given site. For our choice of pa- 
rameters, A = 5 happens to fall in the region of £ > r maxi 
while for A = 2.5, £ < r max . In the latter case we clearly 
see (FigJ^(b), A = 2.5) the (statistical) repetition of the 
network beyond the length scale £. 



% = 25 



1 = 5 




FIG. 1. Spatial structure of connectivity network, (a) 
shown is the typical map of links for a system of 50 x 50 
sites generated from connectivity distributions with A = 2.5 
and A = 5. (b) shown are shells of equidistant sites to the cen- 
tral one in a lattice of 300 x 300 sites. Note that for A = 5, 
shells are concentric and continuous fractals; but for A = 2.5, 
shells are broken. 



K in the distribution results from the finite number of 
sites in the system. The scaling in the inset in F ig. ||(b) 
confirms the known relation K ~ mR d ^ x ^ 1 ^ ||13.[l4|. 



(a) 



(b) 




FIG. 2. (a) The resulting connectivity distribution ob- 
tained from simulations performed on two dimensional sys- 
tems of size R = 400, A = 2.5 and for several values of A: 
(circles) A=2, (squares) A=3 and (diamonds) A=4; they all 
end at a cutoff k c (A). For this case r max > £• In the inset we 
show scaling collapse using same data. The threshold takes 
place at k c ~ A d j^ and confirms the validity of our theo- 
retical estimations, (b) Power law distribution of site connec- 
tivity in the network is showed for R — 100, A — 10 and for 
different values of A: A = 2.5 (circles), 3.0 (squares), and 5.0 
(diamonds). Note that in all cases the distribution achieves 
its (natural) cutoff K. In the inset we show the corresponding 

d 

collapse supporting K ~ R x ~ 1 . For this case, r max < £. 



The degree distribution resulting from our embedding 
method is illustrated in Fig. |[ In Fig. ||(a), £ < r max and 



The different regimes are summarized in Fig. ^. 
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FIG. 3. This diagram shows the six regions where different 
behavior of the network is found: for region A: Tmax R <C 
B: r max < £ < J?, C: £ < r ma:E < _R, D: £ < i? < r MMl E: 
R < £ < r m ax, F: R < r max < £. The diagram can be 
mapped into only four regions where the cutoff fc c and where 
size effect K are expected. A and B: no cutoff and no size 
effect; C and D: cutoff and no size effect; E: cutoff and size ef- 
fect; F: no cutoff but size effect. The two symbols indicate the 
parameters corresponding to Fig. [lj}, ( full diamond) A = 2.5 
and (full circle) A = 5. 

We now address the geometrical properties of the net- 
works, arising from their embedding in Euclidean space. 
To this aim, it is useful to consider the spatial arrange- 
ment of the networks as measured both in an Euclidean 
metric and in chemical space. The chemical distance I 
between any two sites is the length of the minimal path 
between them (minimal number of links). Thus if the 
distance between the two sites is r, then I ~ r dmin de- 
fines the minimal length exponent d m i n . We will see that 
dmin < 1 (for d > 1), contrary to all naturally occur- 
ring fractals and disordered media. Sites at chemical dis- 
tance I from a given site constitute its l-th chemical shell. 
The number of (connected) sites within radius r scales as 
m(r) ~ r df , defining the fractal dimension df. Likewise, 
the number of (connected) sites within chemical radius 
I scales as m(l) ~ l dl , which defines the fractal dimen- 
sion di in chemical space. The two fractal dimension are 
related: d m i n = di/dj [|0-12|. 

To study df, we compute the perimeter S(r), the num- 
ber of sites that connect the interior cluster of a region of 
radius r to sites outside. The fractal dimension then fol- 
lows from the scaling relation S(r) ~ r * _1 . We focus on 
the regime £ > r max . Consider a shell dr', of radius r' . A 
site of connectivity k' within the shell is connected to the 
outside (to a distance larger than r — r') with probability 



P( k ' > (^) d ), eq. (3). Thus, 

S(r) = / dr'r ld ~ x P(k' > (- 



c(X)A. 



,d-l 



r < A, 
r> A, 



where c(A) ~ 1 + — 1) + 1]. In other words, the 

network is compact, df = d at large distances r > A, 
and super-compact, df = d + 1, at r < A. Results for 
df are presented in Fig. |] and are in good agreement 
with Eq. (7j). The slight slope observed for r > A is 
due to analytical corrections, of order r _1 , to the scaling 
S(r) ~ r d_1 , and can be obtained from a more careful 
analysis of Eq. (Q) . 




FIG. 4. Plot of scaled perimeter as a function of the Eu- 
clidean distance from the central site, for several values of A: 
A = 3.0 (top), 3.25, 3.5, 3.75 and 4.0 (bottom). The simu- 
lations where performed with A = 7. Note that the position 
where the curves split, r ~ A, is consistent with our analyti- 
cal results (Eq. (0)). Also, the asymptotic values shown for 
large r are consistent with c(A)A 

In order to compute d m i n (or d{), we regard the chemi- 
cal shells as being roughly smooth, at least in the regime 
i > Tmax, as suggested by Fig. lb (A — 5). Let the width 
of shell I be Ar(Z), then 



I 



dl 



dr 
Ar(l) 



(8) 



since Al = 1. The number of sites in shell I, N(l), 
is, on the one hand, N(l) ~ r(^) d ~ 1 Ar(l). On the 
other hand, since the maximal connectivity in shell / 
is K(l) - W(0 1/(A_1) , the thickness of shell (I + 1) 
is Ar(l + 1) which is determined by the length of the 
largest link to the next shell i.e., r[K(l)], and thus, 
Ar(l + 1) ~ r\K(l)] ~ AK(iy/ d . Assuming (for large 
I) that Ar(l + 1) - Ar(l), we obtain 



.4 



.)«*) 



(7) 



Ar(l) — r^- 1 )- 1 . 
Using this expression in (g), yields 



(9) 



3 



mm ~A-l-lAf {W) 

Thus, above d = 1, the dimensions d m i„ and di = d m i n df 
are anomalous for all values of A. 

In Fig. ||a we plot di as measured from simulations, 
and compared with the analytical result Eq. (|l0|). The 
scaling suggested in Fig. |b, N(l) ~ l d i- 1 $(r r /R d ), is 
valid only for £ > r max . For i? — > oo, we expect that 
the network is scale-free up to length scale £ and the 
analogous scaling will be N(l) ~ Z di_1 \E r (Z di /£ d ), where 
*(i » 1) - a;^-*)/*. 

(a) 



In summary, we propose a method for embedding scale- 
free networks in Euclidean lattices. The method is based 
on a natural principle of minimizing the total length of 
links in the system. This principle enables us to em- 
bed the scale-free in Euclidean space without additional 
external exponents such as assumed by Manna and Sen 
]l6[ and Xulvi-Brunet and Sokolov JlTj . We have shown 
that while the fractal dimension df of the network is the 
same as the Euclidean dimension, the chemical dimen- 
sion di > df for all values of A, yielding d m i n < 1 for all 
A and d > 1. 




(b) 
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FIG. 5. (a) The chemical dimension di as a function of 
A. Note the good agreement between theoretical estimations 
(continuous line) and simulations results (full squares), (b) 
The shape of the &(l d ' /R d ) scaling function is shown for A = 4 
and several lattice sizes: R=1000 (circle), 2000 (square), 2500 
(diamond) and 3000 (triangle ). 
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